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The vector mode of cosmological perturbation theory imprints characteristic signals on the weak 
leasing signals such as curl- and B-modes which are never imprinted by the scalar mode. However, 
the vector mode is neglected in the standard first-order cosmological perturbation theory since 
it only has a decaying mode. This situation changes if the cosmological perturbation theory is 
expanded up to second order. The second-order vector and tensor modes are inevitably induced by 
the product of the first-order scalar modes. We study the effect of the second-order vector mode on 
the weak lensing curl- and B-modes. We find that the curl-mode induced by the second-order vector 
mode is comparable to that induced by the primordial gravitational waves with the tensor-to-scalar 
ratio r = 0.1 at £ « 200. In this case, the curl-mode induced by the second-order vector mode 
dominates at £ > 200. Furthermore, the B-mode cosmic shear induced by the second-order vector 
mode dominates on almost all scales. However, we find that the observational signatures of the 
second-order vector and tensor modes cannot exceed the expected noise of ongoing and upcoming 
weak lensing measurements. We conclude that the curl- and B-modes induced by the second-order 
vector and tensor modes are unlikely to be detected in future experiments. 


I. INTRODUCTION 

The recent remarkable developments of cosmological observations such as the cosmic microwave background (CMB) 
or large-scale structure help us to build the standard ACDM cosmology. The new era of high precision cosmology 
makes it possible to acquire rich information about the expansion history of the Universe or the features of density 
perturbations CHi. It is very important to combine several types of observations to reduce degeneracies between 
cosmological parameters. The weak lensing effect is a key observable for revealing the late-time evolution of density 
perturbations. 

The weak lensing effect can be roughly classified into two observables (for reviews, see e.g., diS])- One is called 
CMB lensing, which is the gravitational deflection by the foreground large-scale structure. In CMB experiments, we 
can measure the deflection angle of CMB photons from observed CMB maps through the reconstruction technique [3- 
0. The CMB lensing signals have been precisely detected by the Planck satellite [ll| and are available to constrain 
cosmological parameters. Next-generation CMB observations are planned [H, [i3|: and CMB lensing will become a 
more important observable in the near future 0. The other observable is called the cosmic shear, which can be 
measured by observing deformed galaxy images. The photons emitted from galaxies are deflected by forward density 
perturbations, deforming the intrinsic shape of galaxies. Ongoing and upcoming imaging surveys such as the Dark 
Energy Survey (DES) [l^, Subaru Hyper Suprime-Cam (HSC) 0, Square Kilometre Array (SKA) 0, and Large 
Synoptic Survey Telescope (LSST) [0, can provide us with high-precision cosmic shear data. Thus, the weak lensing 
survey is becoming a more interesting and active area of measurement. 

The first-order cosmological perturbation theory includes three independent modes: scalar, vector, and tensor. 
Among them, the scalar mode is the dominant component in our Universe and has been well determined by cos¬ 
mological observations. Conversely, the vector and tensor modes are subdominant and have not been observed by 
current observations. In particular, the vector mode is often treated as the negligible component since it rapidly 
decays in the standard first-order cosmological perturbation theory with perfect fluids. Nearly all inflation models 
predict primordial gravitational waves (PGW). With a nonvanishing amplitude, namely a nonzero tensor-to-scalar 
ratio r, primordial gravitational waves correspond to the tensor mode. On the basis of current observations, pri¬ 
mordial gravitational waves have the small tensor-to-scalar ratio of r < 0.1 [3,0. In the context of scalar, vector, 
and tensor decompositions, the weak lensing effect also can be associated with each mode. The deflection angle for 
CMB lensing can be written in the gradient of the scalar lensing potential (gradient-mode) and the rotation of the 
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pseudoscalar lensing potential (curl-mode). The deformation of the shape of galaxies is described by the Jacobi map, 
which can be decomposed into even and odd-parity modes (E- and B-modes, respectively). The vector and tensor 
modes, rather than the scalar mode, induce the curl- and B-modes. Therefore, the weak lensing curl- and B-modes 
are key observables for exploring subdominant modes. 

Some possible sources for the vector and tensor modes in extensions of the standard ACDM cosmology are available. 
The weak lensing induced by the primordial gravitational waves has been well studied [io, [Hi. The primordial 
gravitational waves with r = 0(0.1) do not have detectable amplitudes for the curl- and B-modes, even under the 
assumption of ideal experiments. Cosmic defects are also possible sources of the vector and tensor modes. The weak 
lensing effect induced by cosmic strings has been studied, and weak lensing measurements can constrain parameters 
related to cosmic defects [H, However, in the first-order cosmological perturbation theory, the vector and tensor 
modes must have model parameters, e.g., the tensor-to-scalar ratio or the strength of the cosmic string tension. The 
amplitudes of the weak lensing signal induced by the above models depend on the model parameters and the generated 
weak lensing signal has uncertainties. 

In the second-order cosmological perturbation theory, the second-order vector and tensor modes are naturally 
induced by the product of the first-order scalar modes. These modes do not have free parameters, since the first-order 
scalar mode is well determined by current observations. The secondary CMB polarization anisotropy induced by these 
modes has been discussed in the literature (23 - [HI| . The application of these modes to the weak lensing is also possible 
and quite interesting. The contributions of the second-order vector and tensor modes to the gradient- and E-modes are 
investigated in [2^32[. As the first-order scalar mode can induce the gradient- and E-modes, the contribution from 
the second-order vector and tensor modes to the gradient- and E-modes must be smaller than that from the first-order 
scalar mode. In Ref. [s^, the authors estimated the curl- and B-mode signals induced by the second-order tensor 
mode for the first time. The effect of the second-order tensor mode on the B-mode signal is comparable with that of 
the primordial gravitational waves with r = 0.4 and dominates on small scales, 10 <£. However, the second-order 
tensor mode tends to have a smaller contribution than the second-order vector mode l29l - l3^ . Therefore, the weak 
lensing signal from the second-order vector mode is expected to exceed that from the second-order tensor mode. 

In this paper, we focus on the weak lensing curl- and B-modes induced by the second-order vector mode. The weak 
lensing curl- and B-modes are generated not by the scalar mode but by the vector and tensor modes. Therefore, the 
curl- and B-modes are good tracers of the subdominant mode in the current Universe. As the second-order vector mode 
must have a larger amplitude than the second-order tensor mode, it is important to estimate the weak lensing signal 
induced by the vector mode. This paper is organized as follows. In Sec. H, we review the second-order cosmological 
perturbation theory limited to the vector mode. The standard cosmology comprises the Einstein-Boltzmann system. 
We expand the Einstein-Boltzmann system up to the second order and show the evolution and spectrum for the second- 
order vector metric perturbation. Furthermore, we discuss the details of the second-order vector metric perturbation 
and reveal analytic explanations. In Sec. HI, we summarize the full-sky formalism of CMB lensing and cosmic shear. 
In Sec. IV, we present our main results on the weak lensing signal and offer some discussion concerning detectability. 
We assume four survey designs for the cosmic shear measurement. Finally, we present conclusions in Sec. V. 

Throughout this paper, we use the units in which c = fi = 1 and a metric signature of (—, -f, -f, -f). We obey the 
rule that subscripts and superscripts of greek and latin characters run from 0 to 3 and from 1 to 3, respectively. 


II. SECOND-ORDER PERTURBATION THEORY 

The second-order cosmological perturbation theory is well established in the context of the secondary CMB 
anisotr opy 1 ^ [HI [S^Hll- In this section, we formulate the second-order cosmological perturbation theory, fol¬ 
lowing |43l. l44| . We work in the Poisson gauge perturbed from the flat Friedmann-Lemaitre-Robertson-Walker metric, 
whose line element is written as 


ds^ = a^(?7) -f , (1) 

where r] is the conformal time and 0 ( 77 ) denotes the conventional scale factor of a homogeneous and isotropic universe. 
Moreover we have introduced and Sg^i, as the conformal flat four-dimensional metric and the small metric 
perturbations, respectively. Throughout this paper, we adopt the line element in a spherical coordinate system 
as 


gfj,^dx>^dx‘ 


dij'^ + dx" + X^^abde^^de^ , 


( 2 ) 
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where x is the comoving distance, and Wabdd^dd^ = d6^ + sin^ dd(p'^ is the metric on the unit sphere. The perturbed 
metric in the Poisson gauge can be given by 

<5^00 = -2^ , (3) 

Sgoi = cTi , (4) 

^9ij ~ 2$dij hij . (5) 

Under the Poisson gauge, the vector mode <Ti and the tensor mode hij obey the transverse condition i = 0 

and traceless condition = 0 due to the gauge condition. Here we denote a comma as a spatial derivative and 
raising or lowering indices is done by Sij. 

The first-order vector mode is usually neglected in the standard ACDM cosmology with perfect fluids since the 
vector mode is only decaying in the first-order cosmological perturbation theory. Furthermore, the first-order tensor 
mode would have a small amplitude, namely, tensor-to-scalar ratio r < 0.12 M- Throughout this paper, we therefore 
neglect the first-order vector and tensor modes. Note that, in the first-order cosmological perturbation theory, the 
scalar, vector, and tensor modes evolve independently and it is sufficient to solve the scalar mode only in the first- 
order theory. On the other hand, in the second-order cosmological perturbation theory, the scalar, vector, and tensor 
modes are no longer independent modes. For instance, the second-order vector and tensor modes are excited from 
the product of the first-order scalar perturbations even if the first-order vector and tensor modes are absent. We 
expand the metric perturbation as 'h = -I- -|- ai = and hij = ^h[^\ where we 

have neglected the first-order vector and tensor modes. In the following subsections, we will discuss the standard 
cosmological perturbation theory, which contains the Boltzmann and Einstein equations. 


A. Boltzmann equation 


We consider the Boltzmann equation for the distribution function, where we denote the canonical 

momentum dx^/dA by and A is the affine parameter. In this subsection we present a brief review of the second- 
order Boltzmann equation, following [d^. The distribution function for photons obeys the collisional Boltzmann 
equation since photons interact with electrons via the Compton scattering. The collisional Boltzmann equation for 
photons is written as 

^(x^P^) = C[/], (6) 

where C[f] is the collision term due to the Compton scattering between photons and electrons in the case of the 
photon distribution function. For massless neutrinos, the collision term in their Boltzmann equation must vanish. In 
order to describe the perturbed Boltzmann equation, it is useful to change the coordinate system from the Poisson 
gauge to the local inertial frame. Hence we use the the momentum in the local inertial frame = {E,ph) such that 
fdhi = 1 rather than that in the Poisson gauge P'^, hereafter (see [H). Here E is the energy measured in the local 
inertial frame, which obeys the Einstein relation E^ — = m^. With these variables, we then expand the photon 

distribution function up to the second order as 

f{r),x,p,n) = f‘'°\r),p) + f^^\r),x,p,n) + ^f^^\r),x,p,n) . (7) 

The zeroth-order distribution function f^^'^{r],p) is fixed to the Planck distribution. It is convenient to define the 
brightness function {rj, k,n) in Eourier space as 


{rj, fc, n) 


J dp p^f<^^’^\T],k,p,n) 
J dp p3/(0)(7y,p) 


( 8 ) 


where we have translated from real space to Eourier space. Furthermore, we expand the brightness function using 
spherical harmonics to eliminate the angular dependence: 




(9) 


The brightness function is related to the temperature perturbation of the photons . Here, the multipole coefficient 
Ai^rn obeys the hierarchical Boltzmann equation 




Ci+l^m A (1,2) Cl^rn « (1,2) 

2£ + 3 ^+1’'" “ 27^ 




( 10 ) 
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where a dot denotes a derivative with respect to the conformal time, = \/P — rn?, and 
multipole coefficients of the collision and gravitational terms, respectively. The multipole coefficient of the collision 
term corresponds to the right-hand-side in Eq. whereas the gravitational term comes from the left-hand-side in 
Eq. ([S]), i.e., the perturbed geodesic equation for photons. These explicit forms are written in [i^. Although in the 
first-order perturbations the scalar (m = 0 ), vector (m = ± 1 ), and tensor (m = ± 2 ) modes are completely decoupled 
each other, the vector and/or tensor modes can be generated due to the nonlinear interactions of the scalar ones when 
the second-order perturbations are taken into account. Note that the multipole coefficient of the neutrino distribution 
function A// ^ can be formulated in the same manner without the collision term. 


B. Evolution equation for second-order vector mode 

We denote this subsection to explicitly write down the evolution equation of the second-order vector mode, which 
can be derived through the space-space components of the Einstein equation. The second-order Einstein tensor, 
and energy-momentum tensors for the relativistic particles T/j and nonrelativistic particles T^j are given by 


and 








-f 


^$(2)^ 4 /( 2 ) terms^ -I- (diagonal part) S^j , 


-I- (diagonal part) <5*^ , 

+ (diagonal part) S^j . 


( 11 ) 


( 12 ) 

(13) 


where p,v, and 11 are the energy density, velocity, and anisotropic stress, respectively. According to the results of 
the linear perturbations theory, we can split the first-order quantities into the primordial fluctuations and 

the transfer functions as $(^^( 77 , fe) = (j)P> {k)^'T{krf)^ ( 77 , fe) = (fc) 4 'T(fc 77 ). The primordial fluctuations in this 

paper are assumed to be the random Gaussian field characterized by the primordial power spectrum. 


(</(i)*(fc)(/(i)(fe')) = {2^fP^{k)5l{k - k') . 

We will further assume a power-law spectrum in the form 


(14) 


27r^ 


pm 



(15) 


where /cq denotes the pivot scale. Hereafter we adopt the scale-invariant spectrum with A'^{ko) = 2.4 x 10 ® [T| and 
TTg = 1 for simplicity. The velocity perturbations can be related to the scalar metric potentials through 


In order to describe the evolution equation for the second-order vector mode, it is convenient to expand the 
variables in terms of the mode functions. To do this, we adopt spherical coordinate basis vectors which are written 
in the Cartesian coordinate system as 


n = (sin 0 cos (/?, sin0 sin (/?, COS0) , (17) 

ee{n) = (cos0cos(^,cos0sin(/7, — sin0) , (18) 

e^{h) = (—sin0sin(^,sin0cos(/5,0) , (19) 

where we have imposed e° = 0, n • = 0, and Ba ■ ef, = uiab in the observer rest frame. By using these basis vectors, 

we define the polarization basis as 



(n) = eg{n) ± 


( 20 ) 
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We then introduce the operator which can be defined in terms of the polarization vectors in Fourier space, 

e(=*=^(fc), as (see, e.g., Refs. 


0^^\k) = , ( 21 ) 

with k = k/k. Since these operators satisfy the transverse condition, i.e., = 0, the vector metric pertur¬ 
bations can be expanded in terms of as 

cr,(? 7 ,a;) = / ax{r],k)0^^\k)e-^'^'^ . (22) 

^ ' A=±l 


With these conventions, we can obtain the equation of motion for the second-order vector mode by projecting the 
Einstein equation as 

df (fc) + 27^af (fc) = 5f (fc) , (23) 

where denotes the second-order source terms defined as 


5f)(fc) = 


+87rGaV(,°^A4,A(^)) 

, f d?ki 

-4/ci 


ci>W(fci)vI/(i)(fc2) . —KM 


(27r)3 

(Pki 4 kl r 
(27r)3 ^3 k I 


1 / ZLtt 




m=b,dm 


dPki 

(^ 


^mo(^l)^mo(^ 2 ) 


( 1 ), 


3-b:o(fci)^f b:A(fc2), 


(24) 


where k = ki + k 2 , for simplicity, we omit the time dependence in the above equations, and we have translated from 
the real space to Fourier space. We can see that if we consider the linear theory with perfect fluids, namely, without 
the right-hand-side in Eq. (l24l) . the vector mode decays with oc a“^. However, in the second-order cosmological 
perturbation theory, the nonlinear mode coupling between the first-order scalar modes induces a nondecaying vector 
mode, which is an entirely second-order effect. Consequently, the second-order vector mode imprints on the weak 
lensing. 

Before showing the numerical results of the second-order vector mode, we note the features of our numerical 
calculation [i^, To solve the second-order Einstein-Boltzmann system Eqs. and (1251) . we need the first- 
order perturbations, e.g., or in the Poisson gauge. We obtain these first-order variables by using the 

public Boltzmann code such as CAMB [d^. Furthermore, we store the first-order variables in fc-space from fcmin = 
5.0 X 10“3 [/iMpc~^] to fcmax = 5.0 X 10^ [/iMpc“^]. In addition, we solve the Boltzmann hierarchical equation (ITUl) 
at ^ = 30 for the first-order Boltzmann equation and at £ = 25 for the second-order Boltzmann equation. We checked 
that the results are stable against these choice. 


C. Power spectrum for second-order vector mode 


Before evaluating the effect on the weak gravitational lensing, in this subsection we show the resultant second-order 
vector mode by performing a fully numerical calculation and we then discuss the feature of the power spectrum under 
some approximations. To do this, let us define the unequal-time power spectrum for the vector mode as 

{al{r],k)ay{r]',k')) = {2TrfSxyS^{k-k')P^{r],r]',k) , (25) 


where (• • •) is the ensemble average. We derive an expression for the power spectrum by solving the evolution equation. 
Equation (1251) is easily integrated as 




Kv) 


drj' 




(26) 
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FIG. 1; The evolution of the second-order vector metric perturbation for scales from fc = 10 ^ hMpc ^ to 10^ /iMpc ^ as 
indicated in the figures. 
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Hence we have 

^V2a^ {vi)a^ {V2) , k)S^y (?? 2 ,fe')) ■ (27) 

Once we obtain the brightness functions for photons and neutrinos by solving the Boltzmann equation (jlOll and 
substitute the first-order results for the scalar metric potentials into Eq. (I24E we can obtain the power spectrum 
for the second-order vector mode though Eq. 1221 ). We now solve the evolution equations for the vector mode by 
performing a fully numerical calculation. Figure [T] shows the equal-time power spectrum for the vector mode induced 
by the second-order source terms. 

For illustrative purposes to show the dependence on the wavenumber, we adopt the wavenumbers as from k = 
10“"^ hMpc”^ to 10^ /iMpc”^. The resultant power spectrum for the second-order vector mode during the radiation- 
dominated era seems to grow as oc a on super-horizon scales, while it decays on small scales after it enters the horizon 
scale. This is because the source of the second-order vector mode, namely the scalar potential, decays during the 
radiation-dominated era on sub-horizon scales. In contrast, during the matter-dominated era it always evolves as 
oc 0 ( 77 ) for those wavenumbers. Therefore the second-order vector modes that enter the horizon after the matter- 
radiation equality time do not undergo the above suppression. In Fig. [2l we plot the dimensionless power spectrum 
with various values of the redshift. We find that it scales as on large scales and k~'^ on small scales and its peak 
would be determined by the time of the matter-radiation equality. We study the analytical description of the power 
spectrum of the second-order vector mode in the next subsection. 


D. Analytic description of the power spectrum 

In this subsection, we investigate the feature of the power spectrum for the second-order vector mode analytically. 
It is difficult to estimate the second-order vector mode analytically including the purely second-order quadrupole 
moments for photons and neutrinos. However, if we assume that the purely second-order quadrupole moments for 
photons and neutrinos, and give negligible contributions, the second-order vector metric perturbations 

are sourced only from the convolution of the first-order scalar metric potentials. Indeed, numerical computations 
reveal that the corrections of the purely second-order quadrupole moments to the vector mode amount to only about 
0 ( 10 “^)% during the matter-dominated era, whereas, the quadrupole moments contribute several tens percent 
to the vector mode during the radiation-dominated era. Even in the radiation-dominated era, the scalar potentials 
are still dominant in Eq. (1241) . In this paper, we focus on the weak lensing signals, which are mainly determined by 
the contributions after the matter-radiation equality, in which regime the quadrupole moments contribute at most 
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k 

oc k^ 

1 + z = 10^ 

1 + 2 = 10 "^ . 

1 + z = 10^ . 

1 + ^ = 10 ^ - 

1 + z= 10^ . 


FIG. 2: The spectra of the second-order vector metric perturbation for redshifts from 1 + z = 10® to 1, as indicated in the 
figure. Before matter-radiation equality, the feature of the second-order vector metric perturbation was determined by the 
horizon scale at each time. On the other hand, after matter-radiation equality, it was determined by the matter-radiation 
equality scale, namely, fceq ~ 10“^ /iMpc“^. We can see that the evolutions are same for all scales after matter-radiation 
equality 1 + Zeq ^ 3.3 x 10®. 


0(10“®)%. Therefore, it is sufficient to consider only the scalar metric potentials and we ignore the second-order 
quadrupole moments, if we give a rough estimation. 

To simplify the analysis, we adopt the condition such that the two first-order scalar metric potentials are equal, i.e., 
While this is valid only if the first-order quadrupole moments are negligibly small, we keep it just for a 
qualitative understanding of the behavior of the power spectrum for the vector mode. During the radiation-dominated 
era, the scalar potentials are constant on super-horizon scales while they decay on sub-horizon scales. On the other 
hand, they freeze on all scales during the matter-dominated era (e.g., see 0). We note that under this condition the 
second term of the right-hand-side in Eq. (l24l) should vanish. This is understood as follows: The condition we impose 
here implies that fei and fe 2 are interchangeable. Moreover, the spherical harmonics has the following property: 





(28) 


where we have imposed k = ki + k 2 . As a result, the second line in Eq. (IMl) gives negligible contributions to the 
vector mode in the absence of the quadrupole moments. 

Let us evaluate the vector mode during the radiation-dominated era. Since the fourth term in Eq. (1241) is estimated 
through Eq. (HU) as STrGa^pm^u^o (Pm^ /p(o))<j) 2 , jg suppressed by the factor 1 compared with the 

third term. Therefore we found that the third term gives a dominant contribution to the second-order vector mode. 
Using the explicit expression for the spherical harmonics, the power spectrum for the vector mode induced by the 
third term in Eq. can be written as 


oc k 


(29) 


where k 2 = k ~ fci, k ■ ki = cosOi, and the integrated transfer function T{ri, fci, ^ 2 ) is defined in terms of the transfer 
functions for the scalar potential as 

T(?7, fci,fc2) = 2 f \ [ d?7'a^(?7')4’T(fci?7')^T(fc2?7') • (30) 

a (??) Jo 

With a help of the definition of k 2 and introducing the direction cosine = cos 9i , Eq. (l29l) can be reduced to 
A® 

oc fc y y d/ri fc®P0(fci)P0(fc2)T®(?7, fci, ^2) (2fci^i - fc) (1 -/x^) 


(31) 






To perform this integration analytically, we assume that the transfer function of the scalar potential during the 
radiation-dominated era is approximated as 


^ (77 < 7?eq) • 

1 -f [kr]) 

Substituting the above transfer function into Eq. O, we have 


(32) 


T(?7, fci,fc2) 


1 

r]'^kik2 {ki 



ki arctan (^ 277 ) ~ arctan {kirj) 


{11 < Veq) ■ 


(33) 


In order to evaluate the behavior of the power spectrum, we split the integral of ki in Eq. (I3ip into two parts: the 
contributions from ki > k and ki < k for given k. In the former case, the dummy variables fci and k 2 are related 
through k 2 = fci [l — (fc/fci) 7 t-|- 0((fc/fci)^)] and the integrated transfer function can be reduced to the following form: 


r(77, fci,fc2) 


1 

2k\rf 


arctan(fci? 7 ) 


kiT] 

1 -f (A:i77)2 


77Tl(fel77) (77 <77eq) . 


(34) 


Hence the contributions from the products of the first-order scalar potentials with their wavelengths shorter than k 
are 

^P<y oc k dki J ^dniklP^{ki)r]^ (Ti(a:i))^ 77 ^( 1 - 77 ?) 

nOO 

ockf] dxi {Ti{xi)f = k-qI3i{kr]) , (35) 

J kr) 

where we have used the scale-invariance of the primordial power spectrum Eq. (ESI and we have changed the variable 
ki to xi = kiq. Since ti{x) behaves as x^ for a; 1 and x~^ for a: ^ 1, the integral of /?i in Eq. (I55|) can be evaluated 
as a function of kq : /?i oc (kq)^ for kq ^ 1, and /3i oc (kq)~^ for kq ^ 1. Substituting this into Eq. (I35|) . we calculate 
the contributions from modes with fci > A: in Eq. dSU): 


—Po.(fc < ki,q < ?7eq) oc 


(kq)^ 

ikq)~^ 


fc77 < 1 
^ 77:^1 


(36) 


We can reproduce the behavior of the power spectrum, namely oc k^ for super-horizon scales and oc k~^ for sub-horizon 
scales, which can be seen in the numerical calculations. 

Eollowing the same manner, we can analyze the opposite case, namely fci < k. Expanding Eq. (1311) in terms of the 
small quantity fci/fc <C 1, we find that the leading order term vanishes due to the angular integration. Furthermore, 
we also find that the-next leading order term is suppressed by the power kf. Hence the contributions from modes with 
their wavelengths longer than k are suppressed by the factor ki/k and can be treated as subdominant components. 
Note that although the above estimations do not work around k ^ ki, the results in Eq. (1551) is expected to be still 
correct as long as we estimate the behavior roughly for the following reasons. The integrand in Eq. m does not 
diverge at k = ki. The estimations on fc < fci and fc > fci are smoothly connected. Therefore, the contribution from 
fc Ri fci is at most same order as that from fc < fci. Combining these results, we conclude that the power spectrum 
for the second-order vector mode during the radiation-dominated era is determined by the convolution of the scalar 
potentials with shorter wavelengths. 

We discuss the peak shift of the second-order vector mode shown in Fig. [2] from 1 -|- 2 ; = 10^ to 10^. During 
the radiation-dominated era, the first-order scalar potential remains constant on super-horizon scales while it decays 
on sub-horizon scales. Therefore, the second-order vector mode can grow due to the constant scalar potential on 
super-horizon scales. On sub-horizon scales, the second-order vector mode conversely decays due to the decaying 
scalar potential. As a result, the peak of the second-order vector mode is determined by the horizon scale at the 
corresponding era and the peak keeps shifting until the matter-radiation equality. 

We next consider the vector mode after the radiation-dominated era. The evolution during this era can be easily 
understood through Eq. (l26l) . Generally, when the second-order source term remains constant (i.e., = const.), 

the vector mode evolves as oc q^. This condition is actually satisfied since the scalar potentials during the matter- 
dominated era freeze on all scales, as mentioned above. Hence the evolution of the second-order vector mode is given 

by 


oc 


oc ocr 7 


(for all scales) . 


(37) 
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During the matter-dominated era, the shape of the spectrum for the vector mode does not dramatically change since 
the growing features are the same over all scales. Therefore, the information about the power spectrum during the 
radiation-dominated era propagates to one during the matter-dominated era, i.e., the dimensionless power spectrum 
during the matter dominated era is still in proportion to for super-horizon scales and k~'^ for sub-horizon scales, 
respectively. Although the global feature can be understood as above, in more detail small shifts of the scalar 
potential such as $ —?► 9/10$ during the matter-radiation equality induce an additional small suppression of the 
second-order vector mode, as seen in Fig. [T] After the universe is dominated by the dark energy, the scalar poten¬ 
tials begin to decay for all scales, implying that the second-order vector mode generated by these potentials also decays. 

Before closing this subsection, we introduce the analytic model of the power spectrum, which is originally derived 
in (hereafter referred to as MHM). The explicit form of the approximate solution can be written as 

= Wv{k/k,)F{z)F{z') , (38) 

where Cy ~ 0.026, Wy{x) = (l-|-5x-|-3a;^)“®/^, and fc* = Mpc“^ with Dmo and h being the present cosmological 

parameter of the nonrelativistic matter and the Hubble constant Pq in unit of 100 [kms“^ Mpc“^], respectively. The 
function of a redshift F{z) is given by 


2g\z)E{z)f{n^{z)) 


(39) 


where E{z) = Di„o(l + z)'^ + (1 ~ f^mo), Dm( 2 :) = Dmo(l + z)^/E'^{z). We adopt flmo = 0.27 as the fiducial value. We 
denote f{D,miz)) and g{z) as the dimensionless linear growth rate and the growth suppression factor, respectively. 
One can find that / and g are well approximated as / ~ flm'^{z) and 


g{z) oc Om(z) 0^/’’(z) - ^a{z) + (1 + Dm(z)/2) (1 -k Da(2)/70) 


(40) 


where Da( 2) = (1 — Dmo) lE'^{z) and we will normalize g so that g(0) = 1 [H, [s^, [ll|. We find that the transfer 
functions derived in MHM and those determined by the numerical calculation match after the matter-radiation 
equality. However, we should emphasize that for the MHM approximate power spectrum the effect from the evolution 
of the vector mode over all wave numbers during the radiation-dominated era is assumed to be neglected. As we will 
see in the subsequent analysis, this approximation leads to the non-negligible difference between the full-numerical 
and analytic power spectrum. 


E. Tensor mode 

As we mentioned in the Introduction, the curl mode of the CMB lensing and the B-mode shear can be generated by 
not the scalar metric perturbations but the vector and/or tensor metric perturbations. In this subsection, to compare 
with the second-order vector mode, tensor modes are considered as alternative sources of the observables we focus 

on. In particular, we consider primordial gravitational waves and second-order tensor mode as intriguing examples 

(± 2 ) 

for tensor metric perturbations. To describe the spectrum for the tensor mode, we define the spin-±2 operator OF ’ 
in terms of the polarization vectors as 


(fc) = - (fe)ef ^ (fc) ■ (41) 

Since this operator obviously satisfies the transverse-traceless condition, the second-order tensor metric perturbations 
can be expanded as 

. (42) 

•I ' a=±2 

With these convention, we define the unequal-time power spectrum as 

{Kiv,k)h„,{r]',k')) = {2nf Saa'Su [k - k') ^Ph{r],g',k). 


(43) 






10 


Primordial gravitational waves are generated in the very early Universe and the representative sources for tensor 
mode. Its effect on the CMB lensing and the shear measurement has been discussed in the literature Sim. For the 
evolution of primordial gravitational waves, we introduce the PGW transfer function {kr]), which basically 

describes its sub-horizon evolution. In terms of this, we can write the power spectrum as 

. (44) 

In our analysis we adopt r = 0.1, nt = 0 as the fiducial values, and use — 3ji{kr])/kr] for simplicity. The 

corrections due to the effects during the radiation-dominated era would be small and we neglect this small correction 
throughout this paper. 

Similar to the case of the vector mode discussed in the previous section, the second-order source terms induce the 
tensor metric perturbations, which are expected to be one of the possible sources of the curl mode and B-mode shear 
signals [s^ . The analytic model of the power spectrum for the second-order tensor mode induced by the product of 
the first order scalar metric potentials has been discussed in 1^ . The approximate form of the power spectrum 
derived in is given by 

with Ct ~ 0.062 and IFt(x) = (l -I- 7x -I- 5x^) The transfer function for the second-order tensor mode is 

V i>^v) = (1 - '^)3oo- 

The correction factor attributed to the effect of dark energy is defined as goo = 9{z) ~ 1.3136. We note that 

this formula is valid only after matter-radiation equality time. While the correction during the radiation-dominated 
era for the second-order tensor mode has been considered in (d^ . for the purposes of comparing the second-order 
vector and tensor modes, we will neglect such correction since it must be small [s^. 


F. Comparison of each model 

For comparison, the power spectra for these models at the present time are shown in Fig. [3l The power spectrum 
for the second-order vector mode is larger than that for the second-order tensor mode. Therefore, it is expected 
that the second-order vector mode induces a larger lensing signal than the second-order tensor mode on small scales. 
Furthermore, on small scales, the second-order vector mode has a larger amplitude than the primordial gravitational 
waves with r = 0.1. In other words, the second-order vector mode has the possibility of being detected by cosmological 
observations on small scales, unlike the primordial gravitational waves. We explain the reason why the discrepancy 
between the exact vector model and the analytical vector model appears. From Fig. [Sj we can see that the amplitude 
of the exact model is smaller than that of the analytical vector model for k > fceq, where fceq ~ 10“^ hMpc“^ is the 
horizon scale at the time of matter-radiation equality. For k > /ceq, the analytical vector model does not consider 
the effect of the small suppression around the matter-radiation equality time. The small suppression on small scales 
arises because the first order scalar potentials, i.e., the source of the second order vector metric perturbations, decay 
so rapidly in the radiation dominated era that the source on those scales can sustain the vector perturbations and 
make them grow proportional to the scale factor only after some time has passed since the matter-radiation equality 
time (see the blue and magenta lines in Fig. 1). This small suppression is not included in the analytical model. The 
exact model is about ten times smaller than the analytical vector model for k > fceq. 

From the above discussions, we can understand when the suppression is determined. The analytical vector model 
does not change its peak since this model is calculated in the fiat ACDM model with matters and the cosmological 
constant. Therefore, the difference between the analytical vector model and exact one is appeared until the matter- 
radiation equality. In conclusion, the factor about 10 suppression is determined at the matter-radiation equality 
time. 

The tendency for the second-order vector mode is quite similar to that for the second-order tensor mode in Ref. . 
In the next section, we will show the numerical results of the weak lensing induced by the second-order vector 
perturbation. 
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k [hMpc^] 


FIG. 3: The power spectrum of the primordial gravitational waves with r = 0.1 (PGW), the second-order tensor mode (Analytic 
tensor), the analytical approximate solution of the second-order vector mode (Analytic vector), and the numerical solution of 
the second-order vector mode (Exact vector) at the present time (1 + z = 1). The second-order vector mode dominates on small 
scales rather than the second-order tensor mode. The second-order vector mode derived by numerical calculation is slightly 
smaller than that derived by analytic approximation on smaller scales. 


III. WEAK TENSING SIGNALS 

In this section, we present a short review of the full-sky formalism for the weak leasing induced by the vector and 
tensor modes following Refs. The weak leasing can roughly be classified into two observed objects. 

First, the CMB photons emitted from the last scattering surface are deflected by the gravitational potentials related 
to the large-scale structure, which is called CMB lensing. The CMB leasing is mainly caused by the scalar gravitational 
potential. However, vector and tensor perturbations can also affect the deflection angle of photons, and the vector 
and tensor modes imprint characteristic deflection patterns on the CMB lensing. The deflection angle of the CMB 
photons can be written as the gradient of the scalar potential (gradient-mode) and the rotation of the pseudo-scalar 
potential (curl-mode). By using the nature of the parity these lensing potentials can be reconstructed independently, 
even when the gradient mode dominates the CMB lensing signals [521 ]. 

Second, photons emitted from galaxies are lensed by the large-scale structure, causing the shapes of galaxies to be 
deformed. This is known as cosmic shear. By studying the deformation pattern statistically, we can distinguish traces 
of the scalar, vector, and tensor perturbations. The deformation pattern of the shapes of galaxies can be decomposed 
into parity-even (E-mode) and parity-odd (B-mode) components. 

In the following subsections, we present the full-sky formalism for the deflection angle and the deformation pattern, 
which are related to the geodesic equation and the Jacobi map, respectively. Note that it is sufficient to work without 
the Hubble expansion since the geodesic equation is invariant under the conformal transformation. 


A. Curl mode 

The projected deflection angle on the celestial sphere in the direction n, Ao(n), is generally decomposed into the 
(parity-even) gradient and (parity-odd) curl modes, expressed as 

Aa(n) = (j){n),a + wln),b e^a , (47) 

where a colon denotes a covariant derivative on the unit sphere and e^a is the covariant two-dimensional Levi-Civita 
tensor. The quantities (j) and zu denote the gradient- and curl-modes of the deflection angle, respectively. Raising 
or lowering indices of the two-sphere vector are done by uiab defined in Eq. (I2|). In our analysis, we focus only on 
the curl-mode as the signal of the CMB lensing. We will expand the curl-mode by the spherical harmonics since we 
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observe the scalar and pseudo-scalar lensing potentials on the celestial sphere: 

w{n) = E (48) 

With the coefficients of the harmonics, the angular power spectrum for the curl-mode is defined by 




1 

2i+l 


i 

'y ^ ■ 

m——l 


(49) 


The deflection angle of the light path obeys the null geodesic equation in the perturbed Universe. Solving the spatial 
parts of geodesic equation for the photon rays at the first-order, we obtain the explicit expression of the curl-mode in 
terms of the metric perturbations as [2^ 

w{riy°‘,a = ~ -X:X^):be''a)| , (50) 

where ys is the comoving distance of the last scattering surface from the observer and Qa in the Poisson gauge is 
written only in terms of the projected vector and tensor modes: 

^aivo-X,Xn) = ^-cr,{r]o-x,xn) + hij{r]o-x,xn)n^'^el{h) . (51) 

Here we used the Born approximation in Eq. (1501) and the basis vectors in real space, n and ea(n), have been 
defined in Eqs. (HZll-dlll). We can obviously see that the curl-mode is generated by the vector and/or tensor metric 
perturbations, as already mentioned in the Introduction. 

The angular power spectra induced by the vector and tensor modes, respectively, are expressed in terms of the 
unequal-time power spectra as 


Cr'^’(v) = ^ / f dx f dx' x)5^,k^> x!)Pa{m -X,Vo- x', k) , (52) 

^ ^0 Jo Jo 

. . ^ poo pxs fXS 

/ dx dx sl]j,{k,x)sl]g{k,x)Ph{r]o-X,Vo-X,k), (53) 

^ Jo Jo Jo 

where are the weight function for the vector and tensor modes, 

<k^>x) = 

^Siik,x) = 

and Per and Ph in Eqs. (15^ and (1551) are the power spectra of the vector and tensor metric perturbations presented 
in Sec. im respectively. 


' (£- ly.jijkx) 
(^ + 1)! X 




2)! j^fcx ) 
’-2)! kx^ 


(54) 

(55) 


B. B-mode shear 


In this subsection, we summarize the formalism for the cosmic shear described by the Jacobi map, which maps 
the intrinsic light bundle to the observed light bundle. The geodesic deviation equation is required to handle the 
deformation of the light bundle [s^. The deviation vector projected on the celestial sphere obeys 


d^e 

dx^ 


r a /-a 


(56) 


where we define the perturbed symmetric optical tidal matrix T“b, which is related to the deformation of the light 
bundle. Its explicit form is given by 




dx dy 


( 57 ) 
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We set the initial conditions in the observer’s frame of reference, C“lx=o “ ^ introduce 

the Jacobi map as 

r=2?“hK- (58) 

This matrix characterizes the deformation of light bundle. Substituting this relation into Eq. (15611 . the Jacobi map 
satisfies the following equation: 


dx^ 




(59) 


where the initial conditions for the Jacobi map are ^ d2J“f,/dx|^=o = From here, we decompose 

the Jacobi map into the spin-0 and spin-2 variables as 

qV = T>af)e“^)e(_) , ±22J = ^abe“^)e(±) ’ (50) 

where we have introduced the polarization basis with respect to a two-dimensional vector on the sky as e^±)(n) = 
e)j_)(n)e“(n). Note that the above polarization basis is in terms of spin-1 variables, namely, —)► exp (±za) 

under the rotation around h by an angle a. Furthermore, we define the reduced shear by using the Jacobi map as 


9 = 


+2 


V 


dT) 


9 = 


-2V 

' oV 


(61) 


With these we introduce the E- and B-modes for the reduced shear fields and present the full-sky formalism for weak 
lensing measurements. We first expand the reduced shear, which is a spin-2 variable, by the spin-2 spherical harmonics 
as 


g{n) = '^{Elra+iBtra)+ 2 Ylm{fl) , 

E,m 

(62) 

i,m 

(63) 


where Egm and Bem have electric and magnetic parities, i.e., (—1)^ and (—1)^+^, and are called the E- and B-modes, 
respectively. The angular power spectrum of these modes is defined as 


^xx' 


1 

2 e + i 


I 


(64) 


where X and X' take E or B. Note that the scalar, vector, and tensor modes can induce the E-mode. On the other 
hand, only the vector and tensor modes can induce the B-mode. 

To obtain the expression relevant for the weak lensing measurements, let us expand Eq. (I59|) and solve it order by 
order. Since the tidal matrix vanishes in unperturbed spacetime, the zeroth-order solution of Jacobi map trivially 
reduces to = xs^“&- Substituting this into the first-order geodesic deviation equation, we obtain the symmetric 

trace-free part of the Jacobi map up to the next-leading order [23|: 

dx^|^|T(77o-X,X«):(af,>-^(xf^(a(?70-X,X«):&>)|+^[/l(a6)]o" > (65) 

where hab = [/]q® = f{r]o — XS^XS^) — /(ryojO), and T is written in terms of the metric perturbations as 

T = — (4' -I- 4)) — cTih® -|- —hijn^W . (66) 

Once multiplying the resultant Jacobi map (1651) by and we can obtain the explicit expression for 

the spin--|-2 and —2 reduced shear fields in terms of the metric perturbations. 

In this paper, we focus on the B-mode shear induced by the vector and tensor modes. Note that we consider galaxy 
observations such as imaging surveys, which does not accurately divide the redshift of each galaxy. To discuss the 
weak lensing measurements from imaging surveys, the redshift distribution of background galaxies N{x) should be 
taken into account. 
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The angular power spectrum for the B-mode shear is written in terms of the weight function and the unequal-time 
power spectrum for the vector and tensor modes, 

^oo /*XS /*XS 

^BB.(v) ^ dx' S^\{k,x)S^\{k,^)Pa(r]Q - x,m - x', (67) 

poo /*XS /*XS 

^BB,(t) ^ dx S^\{k,x)SB\ik,x)Phim-X,Vo-X,k), (68) 

^ Jo Jo Jo 

where the weight functions for the B-mode shear, are given by 


<?('") _ * 
'^B,^ — " 


1 /(£ + 2)!(£-l)! 


! /'° 


c(t) _ 

'^B,^ — 


1 l{i + 2)\{£-l)l 
dV (£-2)!(^+l)! 


dxs 


dxs 


N{xs)ji{kx) 




X 


N{xs)jiikx) 

iVg kx^ 


I Nix) 

4 N, 


i;to) + 2^) 


(69) 


(70) 


and Pa and Ph in Eqs. (l67l) and (l68l) are the power spectra of the vector and tensor metric perturbations presented 
in Sec-ini respectively. In our calculation, we assume a distribution of galaxies 7^(x)j which can usually be taken to 
be (see, e.g., Ref. 0) 


^(xs)dxs =Ng- 


2 ( 0 . 64^„,)3 


exp 




0.64z„ 


3/2' 


d2;s , 


(71) 


where is the mean redshift, and the number of galaxies per square arc-minute Ng is defined as 


Ng = / dxiV(x) . (72) 

^0 

In this paper, we focus on four survey designs: DES [l^, HSC [l^, SKA 0, and LSST [I^. The experimental 
specifications of each survey design are summarized in Table. HI 



/sky 

Zm 

Ag[arcmin 

DES 

0.125 

0.5 

12 

HSC 

0.05 

1.0 

35 

SKA 

0.75 

1.6 

10 

LSST 

0.5 

1.5 

100 


TABLE I: The experimental specifications of DES, HSC, SKA, and LSST. It is shown that the sky coverage /sky, the mean 
redshift Zm, and the number of the galaxies per square arc minute Ag. 


IV. WEAK LENSING INDUCED BY SECOND-ORDER VECTOR MODE 

In this section, we show our main results and discuss the size of the effect of second-order vector modes. We 
now calculate the weak lensing signals from the second-order vector mode by performing the numerical calculation 
(hereafter referred to as the exact vector). For comparison, the results for the signals from the primordial gravitational 
waves with r = 0.1, the second-order tensor mode (analytic tensor), the second-order vector mode (analytic vector) 
are also shown. 

First, we show the angular power spectrum of the curl-mode in Fig.|3]for the CMB lensing measurement. The CMB 
lensing reconstruction technique can decompose the lensing potential into the gradient and curl modes. Although the 
gradient mode dominates the lensing signals, owing to this technique, the information about the gradient and curl 
modes can be extracted independently. Even when we neglect the instrumental noise, we need to take into account for 
the reconstruction noise only. The noise estimated by the ideal CMB weak lensing measurement is determined by a 
cosmic-variance limited reconstruction of the curl-mode i0- We found that the curl-mode induced by the primordial 
gravitational waves dominates on large scales, i < 200, while that by the second-order vector mode dominates on 
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FIG. 4: The angular power spectrum of the weak lensing curl-mode. As we expected in section Hi FI the second-order vector 
mode dominates on small scales. Furthermore, the second-order tensor mode becomes a sub-dominant contribution to the weak 
lensing curl-mode. The expected noise from the cosmic variance limit is also shown. 


small scales, i > 200. As seen in Fig. [31 the power spectrum for the vector mode has a peak at the scale corresponding 
to matter-radiation equality. On the other hand, those for the primordial and second-order tensor modes have their 
peaks at the horizon scales. Therefore, the second-order vector mode can affect smaller scales than the primordial 
or the second-order tensor mode does. As expected, the second-order tensor mode gives a subdominant contribution 
to the weak lensing curl-mode. This feature is similar to the CMB polarization anisotropy and the weak lensing 
gradient-mode [32 1. 

However, unfortunately, even if we consider ideal experiments, i.e., only the cosmic-variance limited error, the weak 
lensing curl-mode signals do not exceed the expected noise. Although the curl-mode induced by the second-order 
vector mode dominates the signal of the curl-mode on small scales, it will be difficult to detect the second-order vector 
and tensor weak lensing signals in future experiments. We conclude that the curl-mode induced by the second-order 
modes cannot be detected by any CMB observations in the future because of the cosmic variance limit. On the other 
hand, recently, a new possibility has emerged of detecting the weak lensing signals in 21cm observations [5^ [s^ . 
The angular power spectrum of 21cm fluctuations can be expanded up to £ ^ 10^ since they do not have diffusion 
scales unlike CMB fluctuations. Furthermore, the 21cm fluctuations enable us to observe the fluctuation at different 
frequencies which corresponding to the different distances. Therefore, the signal-to-noise ratio can be substantially 
improved. For example, in Ref. [5^, the observable scalar-to-tensor ratio reaches r k, 10“®. If this sensitivity is 
reached in the future observations, the 21cm curl-mode induced by the second-order vector mode should be detected. 
The 21cm fluctuations would be a good probe of the weak lensing curl-mode. 

Second, we show the angular power spectrum of the B-mode shear with the four representative imaging surveys, 
DES, HSC, SKA, and LSST, in Fig.jSJ Unlike the CMB lensing, the statistical error in the cosmic shear measurements 
is determined by the intrinsic ellipticity of each galaxy. In this paper, we assume that the error mainly originates 
from the intrinsic ellipticity of each galaxy as 


= 


(Tilt) 


(2£ -b l)/sky 3600fVg(180/7r)2 


(73) 


where (7int)^ ^ is the root-mean-square ellipticity of galaxies. In this paper, we set (7int)^ ^ = 0-3 derived in Ref. [^ . 
The error in the cosmic shear measurements is mainly controlled by the sky coverage /sky and the number of the 
galaxies per square arc minute Ng and we show the error expected by four survey designs in Table. [H From Fig. [SI 
we can see that the B-mode induced by the second-order vector mode dominates on all scales except for the largest 
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FIG. 5: A angular power spectra of the weak leasing B-mode assnmed fonr snrvey designs: LSST (top left), SKA (top right), 
HSC (bottom left), and DES (bottom right). The second-order vector mode dominates the expected signals on small scales. 
The black solid line shows the expected statistical error for each snrvey. 


scale. However, as is the case with the CMB leasing curl-mode, the cosmic shear B-mode induced by the second-order 
vector mode does not exceed the expected noise for each survey design. From Eq. (1751) and Fig. [SJ the combined 
survey design parameter appeared in Eq. (1751) . i.e., y //sky x iVg, should be improved about 10^ compared with LSST 
to detect the B-mode signal. Such an ultimate survey is quite unrealistic even in the distant future in contrast with 
the 21cm lensing observations. 

We note that in this paper, we focus on the standard cosmological model which can characterize the primordial 
power spectrum by the primordial amplitude and the spectral index rig in E g. (1141) . However, the non-standard 
model may enhance the primordial power spectrum on smaller scales (e.g., l59|). The second-order signals are 

sensitive to the enhancement on smaller scales since the mode mixing is introduced by the convolution of the small- 
and large-scale fluctuations. The second-order signals would be useful to probe the small-scale physics related to the 
inflation model. 

Let us consider the difference between the weak lensing induced by the second-order vector and tensor modes. The 
equation of motion for the tensor metric perturbation has the form of a wave equation. Therefore, the second-order 
tensor mode induced by the products of the first-order scalar modes cannot be amplified when the source remains 
constant in the matter-dominated era on sub-horizon scales . On the other hand, the evolution of the vector metric 
perturbation is equivalent to that of the vorticity. We can see that the vorticity with the source is well amplihed in 
Eq. (l26l) even in such an era. Therefore, the amplitude of the second-order vector mode is larger than that of the 
second-order tensor mode. 

To conclude this section, we remark on other second-order contributions to the weak lensing curl- and B- modes. 
Du ring photon propagation, there are some corrections to the weak lensing formula induced by the geodesic effect 
[^[^. The geodesic effect would have the possibility to enhance the curl- and B-mode signals. However, this geodesic 
effect is induced not by the vector and tensor modes but by the product of the hrst-order scalar perturbations such 
as the Weyl potential, which we leave for future work. 
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V. SUMMARY 

In this paper, we explored the weak lensing signals induced by the second-order vector perturbation. The weak 
lensing effects are classified into two observables: CMB lensing and cosmic shear. Both the signals of the CMB lensing 
and cosmic shear can be decomposed into two modes by using parity, namely, the gradient- and curl-modes for the 
CMB lensing and the E- and B-modes for the cosmic shear. The curl- and B-modes are only induced by the vector 
and tensor modes. In the standard cosmology, the vector mode is neglected and the source of the curl- and B-modes 
is limited to the case of primordial gravitational waves, which have not been observed yet. However, when we expand 
the cosmological perturbation theory up to the second order, the vector and tensor modes are naturally induced by 
the product of the first-order scalar perturbations. As the first-order scalar perturbation theory is well established 
by a number of observations, the second-order vector and tensor modes do not include free parameters and are well 
determined. 

We presented the effect of the second-order vector mode on the weak lensing for the first time. The weak lensing 
induced by the second-order vector mode dominates on smaller scales rather than the primordial gravitational waves 
with r = O.I and the second-order tensor mode. In particular, the cosmic shear induced by the second-order vector 
mode dominates on almost all scales. This is because the second-order vector mode can be enhanced when the source 
exists in the matter-dominated epoch while the second-order tensor mode remains constant even if the source exists. 
This difference also affects cosmological signatures such as the CMB polarization anisotropy. However, the weak 
lensing signals induced by the second-order vector mode cannot exceed the expected noise estimated by the cosmic- 
variance limit and the shot-noise for the CMB lensing and cosmic shear, respectively. Therefore, unfortunately, it 
seems difficult to detect the CMB curl- and B-modes induced by not only the second-order tensor mode but also 
the vector mode in the ongoing and forthcoming weak lensing observations. However, the 2Icm observations can 
decrease the expected noise and it may be possible that the 21cm lensing observations can be detect the 21cm lensing 
curl-mode. 

Throughout this paper, we assume the standard cosmological model. In other words, the primordial power spectrum 
is characterized by the amplitude and the spectral index. However, non-standard cosmological models can enhance 
the primordial power on much smaller scales. The weak lensing curl- and B-modes would become the good probe to 
search for the small-scale power spectrum and we leave this to future work. 
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